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Abstract 

We consider a system of partial differential equations describing the steady flow of a 
compressible heat conducting Newtonian fluid in a three-dimensional channel with inflow 
and outflow part. We show the existence of a strong solution provided the data are close to a 
constant, but nontrivial flow with sufficiently large dissipation in the energy equation. 

1 Introduction 

We investigate the stationary flow of a heat conducting compressible fluid in a cylindrical domain. 
The fluid is assumed to be Newtonian. Then the flow is described by the stationary Navier- 
Stokes-Fourier (NSF) system: 



pv • Vi> — div S(Vi>) 


+ V7r(p,9)=pf 


in 


n, 


div (pv) = 




in 


n, 


div (pEv) — pf ■ v - 


- div (irv) + div (S(i>)i>) — div q 


in 


n, 


(S(Vv)n) -T k + av- 


r k = b k , k = 1,2 


on 


r, 


n ■ v = d 




on 


r, 


P = Pin 




on 


r 


-q-n + L(9-T) = 


-- 9 


on 


r, 



Mathematics Subject Classification (2000). 76N10, 35Q30 

Keywords. Steady Navier-Stokes-Fourier system; inflow boundary conditions; strong solution; small data; 



1 



where £1 = £l x [0,1] with fi C M 2 smooth, v is the velocity field of the fluid, p is the 
density, 9 is the absolute temperature. We assume that the pressure w(p, 9) is a twice continuously 
differentiable function on IR + x IR + such that 

7r(l,T ) = Po >0, 
<Vr(l,T ) = p 1 >0, (1.2) 

a fl 7T(l, T ) = p 2 >0. 

We consider the Newtonian compressible fluid, i.e. the stress tensor has the form 

S(Vv) = p(vv + V T v - -divul) + Adivul, (1.3) 

where the constant viscosities p and A fulfill p > and A > 0. 

Further, the friction coefficient a > 0. Finally, g = —kV9 is the heat flux. Note that we 
could treat the situation k = n(p,9), p = p(p,9) and A = X(p,9) with suitable assumptions 
on these functions. It would only lead to further complications, therefore we take them rather 
constant not to hide the main ideas by too many technicalities. 

The boundary of the domain is divided in a natural way into the inflow part T in , the outflow 
part Tout and the impermeable wall r . More precisely, we define 

r in = {x er -.d< o}, 

Font = {x G r : d > 0}, 

r = {x g r : d = o}. 

We consider a channel-like flow, i.e. we further assume that Y in = fi x {0}, T out = fi x {/} 
and T = d£l x [0, 1]. We set L = on Y in U Y out and L to be a given positive constant on 
the wall r . We set g = on T but we admit g ^ on T in U T out (see fig(T]). Under such 

r 




Figure 1 : The domain 



choice, on the inflow and outflow part the boundary condition for the temperature reduces to 
the inhomogeneous Neumann condition. From the point of view of modeling these parts of the 
boundary are artificial and we can measure the parameters of the flow, what can be reflected in 



2 



different values of g. The condition on r means that the flux of the temperature through the 
wall is proportional to the difference between the outside temperature T and the temperature of 
the fluid. The level of thermal insulation of the wall is given by L. Existence and uniqueness of 
the special solution 8 under our choice of boundary conditions is discussed in the next chapter. 
The total energy is given by E = \\u\ 2 + e(p, 8), where the internal energy e is given by 

e(p,8) = e n (p,8) + c v 8 (1.4) 

with 

dMPi°) = \(ir{p,0)-8d e ir{p,8)), (1.5) 



i.e. the internal energy fulfills the Maxwell relation. Without loss of generality we set in ( |1.4[ ) 
c v = 1. Equivalently we can consider system ( |1.1[ ) with total energy balance < | 1 . 1 1 ) 3 replaced with 
the internal energy balance 



We also set 



div(pev) — kA8 = S(Vv ) : Vv — 7rdivt>. (1.6) 
d p e n (l,T ) = ei, d e e w (l,T ) = e 2 . (1.7) 



Known results 

The existence theory for the steady compressible Navier-Stokes equations in the framework 
of strong solutions was intensively studied in eighties, see e.g. Il20ll . lETTl and many other papers. 
In the beginning of this century, inspired by the results of P.L. Lions (see IfTBlO . the attention 
turned rather towards the weak solutions, i.e. towards the existence theory without any smallness 
assumption on the data. The best results in this direction so far can be found in [8J (slip boundary 
conditions) and llTOl (space periodic situation). Note, however, that the theory is not able to treat 
the situation with nonzero inflow/outflow conditions. 

Strong solutions with inhomogeneous boundary data, which are in our scope of interest in this 
paper, have been considered for the first time in [|29ll in Hilbert spaces and later in L p framework 
(CD, (El, El). 

In [|23l a constant flow in the direction of the axis of the cylinder is investigated and existence 
of a strong solutions in its vicinity is shown in the barotropic case. 

Concerning the NSF system, the existence of weak solutions was studied much later, see 
El, Ifl8ll . The best result in this direction can be found in (5]|. Note, however, that similarly to 
the barotropic case, all these results do not allow to treat the inflow/outflow problems. 

The theory of strong solutions for the NSF system describing the thermal effects is much less 
developed, to our knowledge the only result in L p framework is the paper by Beirao da Veiga H 
where existence of strong solutions in a vicinity of a constant flow with zero velocity is shown. 

Special solution. Here we are interested in showing existence of strong solutions, i.e. we are 



looking for (v, p, 8) e W p (Cl) x W p l (Vt) x W p (Q), solving system ( 1.1 ) so that the solution is 
close to a given special solution. Hence we arrive at a problem with some smallness of the data. 
Even though smallness obviously facilitates the analysis, the results of this kind are still miss- 
ing in the theory of compressible flow, especially for heat-conducting fluids (see the overview 
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above). This is mostly due to high complexity of the system (1.1 >i_ 3 , in particular due to its 
mixed, elliptic (or parabolic in time dependent case) and hyperbolic character. Investigation of 
stationary solutions close to given laminar flows is of particular importance if we want to inves- 
tigate stability of special solutions, which is quite well developed for weak solutions but so far 
almost not investigated in strong solutions framework (see [|22l for an overview of results). 

The simplest example of a special solution to perturb is a solution with zero velocity and 
constant density and temperature. Such flow has been investigated in with homogeneous 
Dirichlet boundary conditions. From the point of view of applications it is important to investi- 
gate flows with large velocity, which leads in a natural way to inhomogeneous boundary data. In 
a cylindrical domain a natural example of such solution is a constant flow in the direction of the 
axis of the cylinder: 

v= [1,0,0], p = l. (1.8) 

In other words, we can state that such a flow is natural to be investigated as a special solution to 
the system ( |1.1[ ) in our domain ft 

As for the temperature of the basic solution the simplest choice would be to assume it con- 
stant. It seems however quite artificial with the flow (v, p) defined above. It is more interesting 
from physical point of view to consider a flow with nontrivial temperature function 6 which cor- 
responds to the solution (v,p) as a solution to corresponding energy balance equation. On the 
boundary we assume the temperature of the form T = T + 7\ with constant, possibly large T 
and small variation 7\. Substituting (v, p) to the internal energy balance equation we see that 9 
should satisfy the system 



1 + e 2 )d xi 9 - kA6 = in ft 
*#L + L(p-T -T 1 )=g onT 

Now we write 9 — 9 X + 9 2 — T + (9 — T ) + 9%, where 9 solves 



K g- + L(9-T -T 1 )=g onT. (L9) 



—kA9 = c g in ft, 

^\r_=g onT, (1.10) 
j^f n 9 dx = T , 

and 

c 9 = ^fads. (l.ii) 



Iftl 



r 



The existence of a solution 9q G to ( 1.10[ > is a classical elliptic result, we only need to apply 



symmetry to deal with corner singularities. The details are given in Lemma |6j In particular 9 
satisfies 

I|V 2 £o||l p + \\V9 \\ Lp + \\9 - T \\ Lp < C[\\g\\ w i-i /P(T) + ||cJ Xp(r) ] < C\\g\\ w ^ /P(ry (1.12) 

Now we see that 9 1 satisfies 

(1 + e 2 )d Xl 9 l - kA9 1 = -(1 + e 2 )d Xl 9 - c g in ft, 
K f t + L(9 1 -(T 1 +T -9 ))=0 onT. 
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The existence and uniqueness of solution to the above system is discussed in Subsection 2.2 

Our goal is to show the existence of a solution to < | 1 . 1 [ ) close to (v, p, 9) . Hence it is convenient 
to introduce the following quantity to measure the distance of the data from the special solution: 



Do = ||/|U P (n) + ^ 2 i= i\\h -a^H^i-i/p^ + \\d-nW\\ w 2-i /P{r) 



+ 1 1 Pin - illw^cr™) + IMIv^-^r) + ll T i 
We are now in a position to formulate our main result. 



(D 1 n^iiw^^cr)- 



(1.14) 



Theorem 1. Assume that Dq defined in {1.14) is small enough, k is large enough, L is large 



enough on Tq, a large enough on Tj n and p > 3. Then there exists a solution (v, p, 6) to system 
fl-/.i| ) such that 

\\ v ~ ^llw| + IIP - p\\w$ + ||0 - 9\\w$ < E(D ). (1.15) 
This solution is unique in the class of solutions satisfying ( |7.i5| ). 

Notation. We use standard notation for the Sobolev spaces, i.e. Wp(£T). For the spaces 
defined on f2 we will skip the domain; for example we write L 2 instead of L 2 (0). A generic 
constant that is controlled (but not necessarily small) will be denoted by C, whereas E will 
denote a small constant (i.e. which can be arbitrarily small for the data small enough). We will 
also use the space Loo((0, L), L 2 (^o)), for simplicity we denote it by L O0 (L 2 ). 



2 Preliminaries 

2.1 Auxiliary results 

In this section we collect some standard results which we use throughout the paper. We start with 
the Sobolev imbedding theorem, let us state here the particular case which is used throughout the 
paper (the general version with the proof can be found in JT] Theorem 5.4]) 

Lemma 1. Let fl be defined as above and let f G W^iVt) with p > 3. Then f G Loo(^) and 
there exist C = C(p, Q) s.t. 

u<cimkx. 



Next we recall the Poincare inequality. Among its different versions we will need the follow- 
ing two. 

Lemma 2. Let Y\ and T 2 denote two nontrivial parts of the boundary T = dQ and u G W 2 1 (Q). 
Then 

lkll! 2 ( ri )<^pOl M llL(r 2 ) + l|V«||! 2(f ,)] (2.1) 

and 

IWIIL(^) < ^[ll^llLcrx) + l|V^|||, 2(n) ], (2.2) 

where Cp = Cp(Q). 
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Lemma 3. (interpolation inequality): 

Let 2 < p < oo. Then for all e > 3C(e, p, fi) swc/j f/zaf V/ G M^, 1 (fi): 



L p <e||V/|U p + C||/|| i2 . (2.3) 

Proof. The interpolation inequality in the Lebesgue spaces and the imbedding W p l ^ L q with 

q = oo for p > 3, q < oo, arbitrary for p = 3 and q = for p < 3 yield 



* < c( P ) \\f\\* \\f\\tf < c(p)(\\f\\ Lp + \\vf\\ Lp ) e 11 -'' l|J - h 



L 2 



Now application of the Young inequality leads to (2.3 ). □ 

Another result we use is the following version of the Korn inequality 
Lemma 4. Let Q G C ' 1 , /x > an J if a = f/zen f2 w not axially symmetric. Then 

I fi\Vu + V T u-^divul\ 2 dx + [ (a(u ■ r) 2 + d 2 ) dS > C \\u\\ 2 wl . (2.4) 

Proof. It is a straightforward generalization of [8, Lemma 4]. □ 

We also apply the following well-known result in finite dimensional Hilbert spaces (the proof 
can be found in [28J): 

Lemma 5. Let X be a finite dimensional Hilbert space and let P: X — >■ X be a continuous 
operator satisfying 

3M>0: (P(£U)>0 for U\\=M. (2.5) 
Then 3£* : ||f*|| < M and P(£*) = 0. 

Our proofs will be based on the L p regularity of elliptic problems with Neumann and slip bound- 
ary conditions. These are classical results, in case of cylindrical domain we only have to deal 
with the singularities of the boundary, for the sake of completeness we give the proofs. The first 
result concerns the Neumann problem. 

Lemma 6. Let 1 < p < oo, f G L P (Q), d G Wp~ 1/p (T) and j Q fdx = J r ddS. Let C G R. 
Then there exists a solution u G W 2 (Q) to the problem 

Au = f in Q, 

t= d ° n r > (2-6) 

j^f u udx = C , 

and 

\\u - CbHwj < C [\\f\\ Lp + \\d\\ wr/P{r) ] . (2.7) 
Moreover, this solution is unique in the given regularity class. 
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Proof. In fact, the only difficulty we have to overcome are the singularities of the boundary. If 
d\r in ur out — 0, then we can simply extend the weak solution using symmetric reflection that 
preserves the homogeneous Neumann condition. In the neighbourhood of T in it will be defined 
as 

b »m={ S!: :;<!!; < 2 - 8 > 

where x = (xi, x 2 , x 3 ) and x = (— x 1; x 2 ,x 3 ). 

For Eg (u) we obtain a Neumann problem in a domain with smooth boundary that we solve 
using classical elliptic theory. 

If d does not vanish on T in U T out , we reduce the problem to the previous case. Let us 
consider Y in , on r out we proceed the same way. Introducing an extension u G W^(Q) such that 
7^-|r jn = d and ||mo||w| < C|MllvK 1 " 1/3 '(r)' we cons ider u = u — u that satisfies 



Au = f, 

du I 
dn 1 1 



0. 



where ||/|| ip <C[||/|| ip + ||d||^- 1/P(r) ]. □ 
The second result concerns the Lame system with slip boundary conditions: 

— //div (V« + V t m — div ul) — A Vdiv u = F in f2, 

2//(D(u)n) ■ n + a u ■ n = B u i = 1, 2 on T, (2.9) 

n ■ w = on T. 

The following lemma gives existence of a solution and the maximal regularity elliptic estimate 
that we apply in our proofs. 

Lemma 7. Under the assumptions of Theorem^there exists a unique u G W£(£l) solving (2.9) 
with 

2 

Mw* < C[\\F\\ Lp + HA||^-v P(r) ]. (2.10) 

i=i 



Proof. Under the assumptions on \l and A (2.9) is elliptic so we easily get a weak solution. The 
only problem we encounter showing regularity of the weak solution are the singularities of the 
boundary on the junctions of r with T in and T out . These singularities can again be dealt using 
the symmetry arguments. Hence we define the operator E% s extending a vector field defined for 
{x : X\ > 0} on the whole space as 

where x is as above and u(x) = [— u l (x), u 2 (x), u 3 (x)}. Then we have 



™ v) 



-div {VE: s (v) + V T El s {v) - |div«l) + AVdiv^ sV , 

= -div (Vv + V T v - fdivvl) + AVdivw, y } 
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and on the plane x\ = the extension E v as preserves the slip boundary conditions, i.e. 

2v(D(E v as (v))n) -n + ai E v as (v) ■ n) = 2(i(p(v)n) ■ n + a(v ■ n) 

and n ■ E^ s (v) = n ■ v. Now we can show higher regularity of E^ s (u), extension of the weak 
solution to ( 2.9| ), and since E% s preserves the boundary conditions we conclude that u is a strong 
solution to (2.9) satisfying the boundary conditions on T in . Application of analogous antisym- 
metric extension on T out completes the proof. □ 



2.2 Remarks on the special solution 

Our next aim is to construct solutions to \\.9) . Recall that 

6 = T + (6 -T ) + 6 1 , 



where 9 solves ( | 1 . 1 0[ ) . Therefore by Lemma |6j we have 

re||0o - T \\ W 2 < C(c g + Ibll^i-Vp^) < C|l5 , llw p 1 " 1/p (r)' 
as well as by the Lax-Milgram theorem 

«||#o - T \\ w i < C\\g\\L 2 (r). 



(2.13) 



(2.14) 



Next for 9\ we have problem ( 1.13 ). Our aim is to show existence of a weak solution to this 
problem as well as its regularity. As the problem is linear, the existence and uniqueness of a 
weak solution is direct consequence of corresponding a priori estimates. First we show estimates 
in W\, then in Wp . Since we will require k and L large, we have to control how all constants 
depend on them. For the sake of simplicity we assume that n and L are approximately of the 
same size. 



Multiplying ( 1.13 ) by Q\ and integrating over Q we get 



k J a |Wi| 2 dx + f ro L8 1 dS+l(l + e 2 ) f Vmt 1 dS 
= + e 2 ) J Tin e\ dS + J Fo L9 1 (T 1 + T - 6 Q ) dS - J Q ((1 + e 2 )d xi 9 + c g )Q x dx. 

All terms on the l.h.s. are nonnegative. Applying the Poincare inequality to the first term on the 
r.h.s. and the Young inequality to the other two ones we get 

k /JWi| 2 dx + L J Vq el dS 

< c P ( / ro t x ds + j Q |v^i| 2 dx) + \ f To eUs + i j a \ ve t | 2 dx 



+C(L J \T 1 + T -9 \ 2 dS + f n \V9 1 \ 2 dx + 



Therefore we get, for k and L sufficiently large with respect to C p , 



>' 1 < V9l \ 2 dx + ^ ! el dS<C(\\g\\l (r) + L\\e - ^11^+^11^11^ + ||<9 -T \\l HV) ). 
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Hence, using also (|2.14[) 



PiWwl < C(\\g\\ L2 (r) + ||7i|U 2 (r)) : 



(2.15) 



where C is a fixed constant independent of k and L, provided k ~ L. Uniqueness is a direct 
consequence of ( |2.15[ ); if the data are zero, then the solution is zero as well. 

Finally, proceeding similarly as above we will deduce the existence of strong solutions to- 
gether with the estimates of the corresponding norms. First we rewrite (| 1 . 1 3|) to the form 



kA9 1 = (1 + e 2 )d Xl 9 1 + (1 + e 2 )d Xl 9 + c g in Q, 



.901 
dn 



-L{6 1 + 9 - T - T x 



onT, 



(2.16) 



and using Lemma [6] 



^ II II w2 < C(||^i || w^- + ll^olliVp 1 + IMU 2 (r) 



+ L ll^i|lw„ 1 - 1/P (r) + L H Tl llty„ 1 - 1/p (D + L H^o - T o|| w i-i/p (r) ) 



(r) 



Now, applying the Poincare inequality in the form 

||#i||tyi < £||^i||w2 + C(e)||^i||wi 



and using ( |2.15| ) we end up with 

Lemma 8. Let T x £ Wl~ l ^ p iT), g £ Wp~ l ' p (Y), p > 3 and Ze? k, L be sufficiently large. Then 
there exist a unique weak solution 9\ to (1.13) such that 

Piling < C[||T\|U 2(r) + ||p|U 2(r) ]. 



(2.17) 



Moreover, the solution is strong, i.e. Q\ £ W 7 ^ and 

Pilki < Ctll^llw^-^cr) + I^H^-^cr)]- 
//ere, ?/ze constants C are independent of k and L. 



(2.18) 



2.3 System for perturbations 

As we are interested in small perturbations of (v, p, 9), it is natural to introduce the perturbations 
as unknown functions. For technical reasons it is better to get rid of inhomogeneity on the 
boundary. Hence we construct u £ W% such that u ■ n\ r = d. More precisely, we set u = V(f>, 
where <p solves 

= J T ddS in Q, 
^ = d on r. 

an 

Next we introduce the perturbations 

u = v — v — u , o = p — p, Tj — 9 — 9 — 9 — 9q — 9\. (2.19) 
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Replacing the total energy balance ( | 1 . 1 1 > 3 with the internal energy balance (1.6) we can rewrite 
(fTTT) as 



d Xl u — div S(Vw) + pi Vcr + p 2 V77 




^77) 


in 


Q. 


divw + a Xl + (u + Mo) ■ Va = G(w 






in 




(1 + e 2 )d Xl rj - kAt] + T p 2 divu = 




77) 


in 


fi, 


(S(Vic)n) ■ Tj + cm • r, = z = 


1,2 




on 


r. 


n ■ u = 






on 


r, 


0" = fin 






on 


r 


«£ + ^ = o 






on 


r, 



(2.20) 



where 



F(u, a, 77) =_(1 +_<r)/ - p 2 V(# + 0i) - (d^{a + 1, 9 + 9 X + v) ~ Pi) Va 
(d e ir(a + 1,9q + 9i + ti)- p 2 ) V(ij + 6 + 0i) + div S(Vu ) - d Xl u 
— {u + Mo) ■ V(m + M ) - cr(v + u + u ) ■ V(w + u ), 

G(u, a) = — (er + 1) divu — o"divw, 

H(u,a,rj) = S(V(« + _u )) ^V(u + w ) - T o p 7r(l + a, 77 + # + #i)_divu 
-(1 + a){u + M ) • V(0 O + 9i + V) + Todivu(p 2 - d e n(l + a,9 + 9^ + 77)) 
-(r/ +_(# - T ) + 9 1 )d e n(l + a, 9 _+ 6 1 + r])div (u + u ) 
+d Xl (9 + 9 1 +r ] ) (e 2 - d e e(l + a, 9 + 9 X + 77)) , 

Bi = h- ar\ l) - (S(Vtt )n) • r f - a« • n. 



(2.21) 



The above expressions result directly from substitution of (2.19) to ( |1.1[ ). It is only worth to 
notice that to derive the expression for H we have applied the continuity equation which yields 



div(pev) = pv ■ Ve = pv ■ Vp d p e + pv ■ V9 dee. 



(2.22) 



Now in < |1 . 1 1 ) 3 the first term on the r.h.s. of ( |2.22[ ) can be combined with the 7rdiv v, what finally 
gives the above form of H(-). 

Let us state at this early stage the following estimates on the quantities F, G, H defined in 
flZ2ll : 

\\F(u, <r,ri)\\ Lp < C[(\\u\\ W 2 + \\a\\ K + \\r)\\ wi ) 2 + (\\u\\ w * + \\a\\ w if + D„(l + A)] , 
\\G(u, a)\\ w i < C[(\\u\\ W 2 + \\a\\ w i) 2 + D Q (1 + A 



\H{u,a,T])\\ Lp < C[(\\u\\ W 2 + \\(r\\ w i + \\v\\w*) 2 

^ 3 1 D Q (1 + A 



\ U \\W2 + \\<7\\W$ + \\V\\W^ 



(2.23) 

where A is defined in (1.14). The bounds (2.23) are straightforward, it is enough to observe 
that F, G and H are composed of terms which are either quadratic or cubic with respect to the 
functions u, a, 77 or linear with respect to them multiplied by something which can be bounded 
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by D . Note that in the estimate of terms contained differences of the derivatives of 7r or e we 
use the fact fact that these functions are twice continuously differentiable and hence we may use 
the Taylor expansion theorem. 

We sketch briefly the idea of the proof of Theorem 1 , and hence the structure of the remainder 
of the paper. In Subsection 2.4 we explain how to deal with the linearization of the continuity 



equation, i.e. equation ( |2.20| ) 2 . To show existence of solutions we can apply a method of succes- 
sive approximations and look for the solution as a limit of a sequence 



d Xl u n+1 - S(Vw 



n+l\ 



PlV<T 



n+1 



P2V77 



n+l 



F{u 



divw" +1 + + {u n + m ) • Va = G(u n , a n ) 

(1 + e 2 )d Xl r] n+1 - KAr] n+1 + T p 2 divu n+1 = H(u n , a n , rf 1 ) 



(S(Vu 
n ■ u n+1 

n+l 



n+V 



Ti + au 



n+l 



■ Ti 



B, 



1,2 



K ~dn~ + Lr l 



n+l 



in 
in 
in 
on 
on 
on 
on 



Q. 
Q. 
Q. 

r. 

r, 

ini 

r. 



(2.24) 



In Section 3 we deal with the linear system (3.1 ) showing the a priori estimates. The first step is 



the energy estimate which is derived in quite a standard way. The energy estimate is then applied 
to show the maximal regularity L p estimates, hence estimates in the space where we look for 
the solution (Lemma fT2l). It can be considered as the main result of Section 3. In the end of 



Section 3 we solve the linear system, hence we show that the sequence (2.24) is well defined. In 
Section 4 we deal with the sequence of approximations. Using the estimates for the linear system 



combined with (2.21 ) we show it is bounded in x x and has a contraction property 
in a weaker space W 2 L x L oc (L 2 ) x W 2 . An analogous estimate gives uniqueness in the class of 



solutions satisfying (1.15). In other words we apply a generalization of the Banach fixed point 
theorem to show the convergence of the sequence of approximations to the solution of ( |1.1[ ). 



2.4 Steady transport equation 

In this section we show the well posedness of the operator S : L 2 — y Loo(L 2 ) defined as 



w 



S(h) 



d Xl w + U -Vw = h in V(Q), 
w = w in on r, ; „, 



(2.25) 



where w is a weak solution to ( |2.25[ ) in the sense as introduced in ( |4.3[ ). The solvability of the 
above problem, to which we refer as a steady transport equation, will be needed to show the 
energy estimate for the density in Section 3 and to solve the continuity equation in Section 4. 
The result is given in the following 

Lemma 9. Let h e L 2 {VL) and w in G L 2 (T in ). Let U G (fi) be such that \\U\\ W 2 is small 
enough and U ■ n\r = 0. Then there exists a unique solution w = S(h) G Loo(L 2 ) and 



\\S{h)\\ Lao (L 2 ) < C [|Kn|k(r in ) + INk(fi)]- 



(2.26) 
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The proof follows H23H . For the sake of completeness we recall here the most important steps, 
referring to Il23ll for the missing details. The idea is to get rid of the term U ■ Va introducing a 
change of variables x = ip(z) satisfying the identity 

d Zl =d Xl + U- V x . (2.27) 

We construct the mapping ip in the following 

Lemma 10. Let ||E/||w 2 be small enough. Then there exists a set O C M? and a diffeomorphism 
x = ip(z) defined on O such that Q = ip{0) and (2.27) holds. Moreover, if z n — > z and 
ip(z n ) — > T then n 1 (^) = 0, where n is the outward normal to O. 

Before we sketch the proof let us make one remark. The last condition states that the first 
component of the normal to ^>~ 1 (r , o) vanishes, but since ij) is defined only on O we formulate 
this condition using the limits. It means simply that the image O = is also a cylinder 

with a fiat wall. 

Sketch of the proof of Lemma 10 The identity (2.27 ) means that ip must satisfy 

A natural condition is that ip(T in ) = Y in . Thus we can search for ip(zi, Z2, z 3 ) = ip Z2jZ3 (zi), 
where for all (z2, 23) such that (z 2 , z 3 , 0) G Y in the function ip Z2 ,z 3 (•) is a solution to a system of 
ODE: 

ds^ Z2 , Z3 = l + U\^ Z2 , Z3 ), d s r z2 , Z3 = U 2 (^ Z2 , Z3 ), d s ^ Z3 = U*(ip Z2 , Z3 ), 
^(0) = (0,Z2,z 3 ). ^ • ; 



Using the regularity and smallness of U we show existence for (2.29), hence we have if) defined 
on some O such that Q = ip(0). The smallness of U implies also that ip(z) = z + "4>e( z ), 
where HV'ellw 1 i s small and we conclude that ip is a diffeomorphism. Let us denote <p — ' l P~ 1 - 
Now it is natural to define the subsets of dO as dO = Oi n U 0Ut U Oq where Oi n = T in , 
Oout = {z : z = lim0(x„), x n ->■ T out } and O = {z : z = lim0(x n ), x n ->■ r }. To show 
that n 1 (z) = for z G O it is enough to observe that 

D<4>(z)([l, 0, 0]) = [1 + U\x), U 2 (x), U 3 (x)], 

where x = ip(z). But for x G T the vector on the r.h.s is tangent to T since U ■ n\ ro = 0. 
We can conclude that on O the image in ip of a straight line {(s, z 2 , z 3 ) : s G (0, 6)} is a curve 
tangent to r , and thus Oq is a sum of such lines and so we have n 1 ^) = 0. The proof of lemma 



10 is completed. □ 



Proof of Lemma^j First we define S(h) for a continuous function h as 

r<t>i(x) 

S(h)(x) = w in (0,(f)2{x),(j) 3 (x)) + / h(i/}(s,^2(x),^ 3 (x)))ds. (2.30) 
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The condition n 1 = on 4>(T ) guarantees that a straight line (s, z%, z 2 ) : s G (0, 6) has a picture 
in ft and thus we integrate along a curve contained in ft. It means that S is well defined for con- 



tinuous functions defined on ft and the construction of ip clearly ensures that S satisfies (2.25 ). 



Next we have to extend S on L 2 (ft). To this end we show the estimate (2.26) for continuous h. 
Let Q Xl denote an si-cut of ft and let x := (x 2 , x 3 ). Then by (|2.30[) we have 



\\S(h) \\l 2(Qxi) = J Qxi [w in (0, + Io Mx) K1>(*> <h(?), ds 



< 2\\wi 



L(r„o + Cfn Xl Io Mx) h*M*> U*)M*))) 



dx 



dsdx < C \\\w 



Jin \\L 2 (T irl ) 



+ \\h\ 



L 2 (n) 



]■ 



The above estimate holds for every x\ G (0, L) what implies (2.26). Now we can define S(v 



for v G L 2 (fi) using a standard density argument, the details are given in (" ll23~T0 . To show the 



uniqueness of the solution we can rewrite the r.h.s of (2.25 ) as 



in 



Q. 



W = W ir . 



l +U l 

on r, 



(2.31) 



where the equation is to be understood in the sense of distribution, and, treating x\ as a "time" 
variable, adapt Di Perna - Lions theory of transport equation ([5]) that implies the uniqueness of 
solution to (2.31 ) in the class L 00 (L 2 ). The proof is thus complete. □ 



3 Linearization and a priori bounds 



In this section we deal with the linear system 



d Xl u — div S(Vm) + piVcr + p 2 Vr] 
div u + d Xx a + U ■ Vct = G 
fod Xl V + ndiv u — nArj = H, 



(S(Vtt)n) • Tj + au ■ Ti 
n ■ u = 



Bi 



1,2 



' dn 



+ Lr] = 



in 
in 
in 
on 
on 
on 
on 



ft. 
ft. 
ft. 

r. 

r, 

ini 

r 



(3.1) 



with given functions {F,G,H,Bi,a in ,U) G [L v x x L p x Wp~ 1/p (T) x W}(T in ) x W J 
such that 



U ■ n = on r , U ■ n > on T 



out J 



U ■ n < on T, 



(3.2) 



Note that this system actually represents ( |2.20| ), where we only redefined a few constants. The 
main difficulty lies in deriving appropriate estimates for ( |3.1[ ). We start with the energy esti- 
mates, then applying the properties of slip boundary conditions, Helmholtz decomposition of the 
velocity and classical elliptic theory we derive the L p estimates. 



Apart from solving the system (3.1), the estimates derived in this section are used later, 



combined with (2.23 ), to show the convergence of the sequence of approximations (2.24) 
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3.1 Energy estimates 



Lemma 11. Let (u, a, rj) solve system \3.1\ with (F, G, H, B i} a in , U) G V* x L 2 x L 2 x L 2 (T) x 
^(r^n) x Wp with \\U\\wjj small enough, U ■ n on Y fulfilling conditions ( 3.2 ) above, where V* 
is the dual space to 



0}. 



V = {ve W 2 \n) :vn\ r 

Assume further that k, L and a fulfill the assumptions of Theorem^ Then the following estimate 
holds: 

IMIw-i + Ikllw**) + WvWw* 

< C[\\F\\ V * + \\G\\ L2 + \\H\\ L2 + H^lli^r) + ||a ifl ||x 3 (r jn )]. 



(3.3) 



(3.4) 



Proof. We multiply (3.1 >x by u and integrate applying the following identity: 

f n — div S(u) v dx = f n S(u) : Vt> dx — f T (S(u)n) ■ v dS. 
With application of the Korn inequality and the boundary condition ( 3.1[ >4 we get 
C|I w IIvp 2 — f n ??div u dx < J n F-udx + J Q Gadx 
+ fj n a 2 div Udx + f J r . n al(l -U-n)dS + E-=i J r B iU ■ r< dS. 
Next we divide < |3 . 1 >3 by r± (recall it is constant), multiply by r\ and integrate. It yields 

±J u \W\ 2 dx + ±J ToV *dS 
- ^ fr in V 2 dS+^ J Fout r] 2 dS + f n r]div udx = f n H V dx, 

hence 



(3.5) 



(3.6) 



K 

n Jn 



\Vr)\ 2 dx + 



r 



2r 1 



T) 2 dS. 



(3.7) 



1 ir, 



7] dS + r]div udx < / Hrj dx + 
Jn Jn 

To get rid of the boundary term on the r.h.s we can apply the Poincare inequality ( |2.1[ ) and rewrite 



( gJJ ) as 
1 



K — 



\Vr]\ 2 dx 



1 



L — 



C, 



i] dS + r]divudx< / Hr/dx, (3.8) 



where Cp is the constant from (2.1 ). For k large enough and L large enough on T the first two 



terms on the l.h.s. will be positive. Now we can combine ( |3.6[ ) and ( |3.7[ ) obtaining 

C7(ft, k, L) [||w||^ + J Q \Vr]\ 2 dx + \\v\\l 2i r )} < J n H V dx + J Q F ■ udx 
+ j n Gcr dx + \ f u a 2 div U dx + \ / r , of n (l — U ■ n) dS + J2f=i Ir ^ M ' r * ^ 
Using (2.2), the imbedding W^ 1 L ro and Holder's inequality we derive 

IMI^ + \\v\\ 2 w i < \\F\\vh\\ W } + \\G\\l 2 \W\\l 2 

L(r in ) + ^ll°"lli 2 + Yh=1 \\ B i\\L2{T)\\u\\L 2 {T)- 

Recall that E denotes a small constant and so to complete the proof we have to find the bound 
on ||o'||i 00 (i 2 ). But this follows from the existence result in the previous subsection, see ( |2.26[ ). 
Hence we have 

(3.10) 

□ 



W\\l x (l 2 ) < C(\\u\\ w i + 



\<j 



in\\L2(Ti n ) 



Combining (3.9) and (3.10) we conclude (3.4). 
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3.2 L p estimates 



The main result of this section is 



Lemma 12. Let (u,a,rj) solve system (3.1) with (F,G, H, B i: a ini U) G L p x Wp 1 x Lp x 

Wp 1 ^ p (r) x WpiTin) x Wp, where \\U\\w^ is small enough and K,L,a satisfy the assump- 
tions of Theorem [7] Then 



<C\\\F 



\G\ 



\ u \\w% + w\\w} + WvWw* 



+ \\H\\h p + lkin||tVi(T in ) + £i=i llAllw/i-i/^r)] : ~ Am- 



(3.11) 



The proof of < |3 . 1 1\ will be performed in several consecutive lemmas. We show the bound on 
||r7|| W 2. Then we proceed with the bound on ||c||wi which is the most demanding. First we show 
bound on the vorticity which, together with Helmholtz decomposition of the velocity, makes 
possible to eliminate the term divw from the continuity equation, leading to < |3 .23 >. Using this 
equation we show the bound on the density. Then ( |3.11[ ) easily follows from the classical elliptic 
estimate on |H|iy 2 - 



We start with the bound on rj; to show it let us rewrite (3.1 (3 7 as 

— kAt] = H — rod Xl i] — ridivu in Q, 



drj i 
dn I 



K 1 



r = --// on I . 

The classical elliptic estimate (0, 0) for the above system yields 



(3.12) 



Applying the trace theorem to the boundary term and then the interpolation inequality (2.3 ) we 
get 

Nlw* < C[\\H\\ Lp + \\v\\w$ + Nlw, 1 ] + e Nlw*- ( 3 - 13 ) 



By the energy estimate (3.4) we conclude 

\\ri\\ W 2 < CD Un + e||w||w2. (3.14) 

Now we proceed towards the bound on ||cr||wi- As mentioned before, we start with the estimate 
on the vorticity. Taking the curl of (|3.1[)i we get 



u ■ Ti + ^ 



—fiAa = curl [F — d xi u 
a-r 2 = (2xi - 

div a = 



in Q, 

on r, 

on r, 

on r, 



(3.15) 



where \i denote the curvatures of the curves generated by tangent vectors 7$. In order to show 
the boundary relations (13.151)2,3 it is enough to differentiate (13. 1L with respect to the tangential 
directions and apply (|3 . 1 p 3 . A rigorous proof is given in |[T9ll or [|23l . 
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The condition div a = in results simply from the fact that a = curl u. We introduce this 
relation as a boundary condition p.l5[ >4, that completes the conditions on the tangential parts of 
the vorticity. Note that the boundary conditions (0.151)2,3 gi ve me tangential parts of the vorticity 
on the boundary the regularity of the velocity itself and the data. We will now use this feature of 
slip boundary conditions to show the higher estimate on the vorticity (see [|X4 | ] . [|T5Tl . |26J). 

For the above system we have (see [30], Theorem 10.4): 



\<A\w} < C \\\F\\ Lp + H^mII^ + \\Bi\\ w i-i/v (r) + \\u\\ w i- 1/P ] . 



i=l 



Applying the interpolation inequality (2.3 ) and then (3.4) we arrive at 

IMIwj < C{e)D Un + e||w||w2 



(3.16) 



(3.17) 



for any e > 0, where Du n is defined in < |3 . 1 1\ . Now we introduce the Helmholtz decomposition 
of the velocity (see e.g. Q) 

u = V<p + A, (3.18) 
r = and div A = 0. We see that the field A satisfies the following system 



where 

an 



curlA = a in Q, 
divA = in Q, 
A ■ n = on T. 



(3.19) 



For this system we have (see [27]): \\A\\ W 2 < C \\a\\ w i, what by (3.16) can be rewritten as 

2 



\A\\ W 2 < C(e) [||F|| Lp + W B i\\wt 1/p (r)\ + e Mw* 



i=l 



for any e > 0. Now we substitute the Helmholtz decomposition to ( |3.iP i- We get 

V[-(A + |/i)A0 +pi<j] — F — X1 A + /J.AA - d xl V<j) - p 2 Vr], 
but A0 = div u. We denote 

— (A + |/i)div u + pi a = K. 
Combining the last equation with ( 3.1[ ) 2 we arrive at 

7a + d Xl a + U ■ Va = K, 



where 7 



K 



K 



A + 3/i 



+ G. 



(3.20) 

(3.21) 
(3.22) 
(3.23) 

(3.24) 



Equation (3.23 ) makes possible to estimate the W^-norm of the density in terms of - norm of 
K, which in turn will be controlled by (3.21 ) using interpolation and the energy estimate (3.6). 
First we estimate in terms of K. The result is stated in the following lemma 
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Lemma 13. Assume that a satisfies equation (3.23) with K G Wl. Then 



Mlwpi < G [\\K\\ w i + ||o"m||w^(r <n )]- 



(3.25) 



Proof. In order to find a bound on |klU P we multiply (3.23) by |cr| p 2 a and integrate over f2. 
Integrating by parts and next using the boundary conditions we get 



-yhWl ~ I L ^ U\a\*dx + \ J Tout \a\ p (l + U- n) dS 
<\\K\\ Lp \W\\ P C + l p J Vin \a m \ p (l-U-n)dS. 
With the smallness of ||£/||vf 2 > the above implies 



(3.26) 



C\\a\\ p r < \\K\ 



\<j 



ip-i 



+ C\\a. 



and so 



kill* < C [\\K\\ Lp + \\a in \\ Lp{ r m )]- 



(3.27) 



Now we estimate the derivatives. In order to find a bound on d Xi a we differentiate (3.23) 
with respect to Note that VV is not defined in general, however, for 



u := [1 + U\U 2 ,U 2 



(3.28) 



and a G W p l we may set u ■ Vd Xi a := d Xi (u ■ Vcr) — d x .u ■ Va G L p . 

Hence we can differentiate (3.23 ) with respect to x i5 multiply by \d Xi a\ p ~ 2 d Xi a and integrate. 



Since d r .u = d T .U, we have 

/ (k^crl^^aVcr) da; 

Next, since w • Vd Xi a G L p , we can write 



<\\VU\\r \\V*\\ P r<C\\U 



I _i v in ii II Veil? 

Jjr> II II f * r, II II !_/») 



/ n fi • |^(x| p 2 d Xt aVd Xi adx = ± J Q u ■ V\d Xi a\ p dx 
= — - f n \d Xi a\ p div u dx + ^ J r \d Xi a\ p u ■ ndS 
= -I L \d x M P divUdx - I J Fm \d Xi a m \ p (1 + U') dS + \ \d Xi a\ p (1 + U 1 ) dS. 

For % = 2, 3 we have a inyX% G L p (Y iri ) and hence the above defines the trace of \d Xi a\ p on T out . 
We arrive at 



tII^I^ - I /„ div U \d Xi a\ p dx + i J r<mt l^o-l' (1 + t/ 1 ) ^ 
< Ik.^lk, ll^alli; 1 + I f \d Xt a in \ p (1 + U 1 ) dS + C \\U\\ wl \\Va\\l p . 



(3.29) 



For i = 2, 3 (3.29) gives straightforward bound on ||9 Xi cr||i p . In order to estimate d Xl cr we also 
differentiate (3.23 ) with respect to x\ and multiply by \d Xl a\ p ~ 2 d Xl a. The difference is that d xi a 



is not given on T in . To overcome this difficulty we can observe that on T in equation (3.23) 
reduces to 

7a in + U 2 d X2 a in + U 3 d X3 a in + [1 + U 1 ] d Xl a = K, 
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what can be rewritten as 



d Xl a 



71 [K - 7Cr in - U T ■ V T cr in ] , 



l + U l 

where the lower index r denotes the tangential component. Thus we have 

\\dxM\L p {r m ) < C [||-^|r in |U„(r in ) + Ikinllwjcr*,)]- 



Using this bound in (3.29 1, i — 1, we arrive at the estimate 

ll^lli < C [\\d Xl K\\ Lp R^llC + \\U\\w> \\V<r\\ P Lp + ll^ll^cr*.) + W^fw^J. (330) 



Applying the trace theorem to the term ||-K"||L„(r jn ) an d then combining (3.29) (for x 2 and x 3 ) 
with p.30[ ) we get 



\\Va\\% < C [\\VKUJ\VaWl-; + \\U\\ w ,\\Va\\l p + \\K\\^ + \\a m \\^ {r J . (3.31) 

The term ||£7||vp 2 || Vcr||^ can be put on the l.h.s. due to the smallness assumption and thus by 
Young's inequality 

\\V<t\\l p <C [\\K\\ w i + |K|| W 2 C r itt) ], (3-32) 
what combined with ( 3.27| ) yields 

\W\\wi < C [\\K\\ W x + |km|ki(r in )]- (3-33) 
The lemma is proved. □ 
Now we estimate K in terms of the data. The result is 



Lemma 14. Let K be defined in (3.24). Then W5 > we have 



\K\\ W i<6\\u\\m+C{S)D. 



lini 



(3.34) 



where Du n is defined in (3.11) 



Proof. 



c . Applying first ( |2.3[ ) and then ( |3.6[ ) we get 

ii _ ^ ^ \\VK\\ Lp + C{5 X ) [\\F\\ L2 + ||G|| £a + \\B\\ L . 2{T) ] . 

ip hnnnrl nn \\X7 fCW r Rv (1 9 1 1 wp havp 



K\\l p < 5 

bound on || VK\\ Lp . By (3.21 ) we 



Hence, it is enough to find the 

\\VK\\ Lp < C \\\F\\ Lp + \\G\\ Wf , u „ HWg , 

where u = V0 + A is the Helmholtz decomposition of the 
trace theorem and ( |2.3[ ) we arrive at 

\\VK\\ Lp <C[\\F\\ Lp + \\G\U 
+8i\\u\\ W 2 + C(5_ 



(3.35) 



IIV77IUJ, 



A \\w2 

velocity. Applying ( |3.20[ ), ( |3.14| ) 



+*ilH| w? + c7( < y 1 )[||F||x : 

Combining ( |3.35| ) and p.36[ ) we get ( |3.34| ). 



H\\& + ll-Sll^-V^r) 

ll^lUa + ll-^IUatT)]- 



o"m||w p 1 (r i „)] 



(3.36) 

□ 
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The bound on a now follows directly. Substituting ( 3.34[ ) to (3.25 ) we get 



W\\m < S\\u\\w2 +C(8)D Un , 



(3.37) 



where D Hn is given by < |3 . 1 l\ . Now the only missing piece to complete the proof of ( |3.11[ ) is the 
bound on ||m||vf 2 - To show it note that the velocity satisfies the Lame system 



u Xl — divS(Vu) = F — piVa — p 2 Vr\ in fi, 
(S(Vu)n) ■ Ti + a u ■ Ti = Bi, i — 1,2 on T, 
n ■ u = on T. 

The classical theory of elliptic equations (B2l|,|l3l) yields 

IMIw* < c [D Un + \\u\\ w i]. 



Applying (2.3) to the term ||iz||wi and then (3.6) we get 

\w2 < C Di in . 



(3.38) 



(3.39) 



Proof of Lemma 12 We combine (3.14), (3.39) and (3.37). In ( 331) we choose for example 
(5=2^ where C is the constant from (3.39). Now (3.11) follows immediately. □ 



4 Solution of the linear system 



In this section we solve the linear system (3.1 ) and hence show that the sequence (2.24) is well 
defined. We start with showing the existence of an appropriately defined weak solution applying 
the Galerkin method modified in a way to deal with the continuity equation. Then we show the 
regularity of the solutions for the given regularity of the data. For simplicity let us denote 



u :- 



\1 + U\U 2 ,U 3 ]. 



(4.1) 



4.1 Weak solution 



In order to define a weak solution to (3. 1 ) we recall the definition of the space V (3.3 ). A natural 



definition of a weak solution to the system (3.1 ) is (u, a, if) G V x L 00 (L 2 ) x W\ such that 
1. the identity 

/ {v ■ d Xx u + S ( Vtt) : Vf — pi<7 div v — pzrjdiv v} dx + / a(u ■ r») (v-Tj)dS 

2 

= [ F-vdx + Y] [ Bi(vTi)dS (4.2) 
Jn i=1 Jr 



is satisfied V v G V; 
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2. (3.1 



>2 is satisfied in the following sense: 



au-Vcpdx— I a(f) div u dx — / (f)(G — div u) dx + / cr in (j)dS, (4.3) 



V G C 1 (fi) : 0|r out = 0, where w is defined in (4.1 ), and 
3. the equation 



r / d^rjw + n 
In 



wdiv u + k 



Vw ■ V?7 + / Lwr] dS 
n Jr 



Hw dx 



(4.4) 



holds for all w G W\. 

We introduce an orthonormal basis of V formed by {oji}^. We consider finite dimensional 

(Xi G R} C V. The sequence of approximations to the velocity 



air- 



spaces: v = 
will be searched for as u 



N 



E 



N ~N, 



ooi. Due to the equation (3.1>2 we have to define the 



approximations to the density in an appropriate way. Namely, we set a = S(G — divu 



where S : L 2 (0) — > L OQ (L 2 ) is defined in (2.25 ), with w in = a in . 

Finally, to define the approximations of the temperature we introduce R : H 1 — > H l as a 
weak solution operator to equation (|3.1I)3 ;7 , i.e. we set 



R{u 



(4.4) holds with r] := rj N V w G W\. 



The well-posedness of the operator R is direct as it is just a solution operator for a standard 
elliptic equation. To show it is well defined it is enough to recall the estimate from the proof of 
< |3.4 ). In the same way we show the estimate for (4.4) which gives existence of i] satisfying ( |4.4 ) 
for given u G W\, hence R is well defined. In particular, we have 



\\R(u)\\ wi <C[\\F\\ L2 



(4.5) 



With the operators S and R well defined we can proceed with the Galerkin method. Taking 



u 



N 



A? 



OJi, V 



u k , k = 1...N and a = a N = S(G - divu N ) in (4.2), r) = r) 



R(un), we arrive at a system of N equations 



,iV 



B N (u\uj k ) 



0, k = 1 ... N, 



(4.6) 



where B : V —> V is defined as 



■ N ,v N ) 



-Pt J n S(G -div £ 



J n {«V + S(^:Vt) w }(fa 



N \ 



div v N dx — p 2 fa R(£, N )div v N dx 
Tj) dS — f n F ■ v N dx. 



B; 



(4.7) 



Now, if u N satisfies ( 4.6|) f or k = 1 . . . N and <r N and r] N , defined as above, then the triple 
(yN a N r jN\ oofiofioo /l/i oh /l/i /ik f™- ^ „,A cl (\rN n^lciw ^ u/i M- — n w» ,„;n ™n 



satisfies (4.2)-(4.4) for (w,0,w) G (V* x C 1 ^)) x W\, 



0. We will call 



such a triple an approximate solution to ( |4 . 2| )— ( |4~4] ) . 



The following lemma gives existence of a solution to system (4.6): 
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Lemma 15. Let F,G e L 2 (Q), a in e L 2 (T in ), Bi e L 2 (T), i = 1,2. Assume that \\U\ \w2 and 
k, L, a fulfill the assumptions of Theorem^ Then there exists u N G V N satisfying {4.6) for 
k = 1 . . . N. Moreover, 

|| u N || w i < C(DATA). (4.8) 

Proof. We will apply a well-known tool of the Galerkin method, Lemma [5J Thus we define 

P N : V N -> V N as 



P N {^ N ) = ^2B N ^ N ,u k )u k for t»eV N . 



(4.9) 



We have to show that (P (£ ), £ ) > on some sphere in V . As B is linear in the second 
variable, we have 

{P(Z N ),£ N ) = B N (£ N ,£ N ) = f S(V£ N ) : V£ N dx 

Jn 



h 



+ / t N d xl £ N dx+ / a{i N ■nydS-px / S(G — div ^ iV ) div dx 



(4.10) 



^2 



-P2 [ R{t N )div v N dS- [ F-^dx-Y] [ Bi(£ N ■ n) dS. 
v J J1 . Jn i=i ^ 



By the Korn inequality we have 



h+h>c\\^r wi 



(4.11) 



for a large enough. To deal with J 4 let us denote for a moment rj N = R(£ N ). Then setting 
w = i] N in (4.4) (with ^ instead of u) we have 



Z 4 = -P 2 / ^div^ = ^ 



d xl r] N V N dx + n / |Vr/T^ + / L(rj ) dS — / tfr^ dx 



and repeating the reasoning from the proof of (3.4) we infer 



h>-\\H N \\ L2 \\^\\ L 2>-\\H 



N 



\H 



N 



+ U N \\wt). 



We have to find a bound on J 3 . Denoting a N = S(G — div ^ iV ) we have 

- f a N dw£ N dx= [ a N (d Xl a N + U-Va N )dx- [ a N Gdx. 
Jn Jn Jn 



(4.12) 



(4.13) 



Using (2.26) we get 



- / rj N G N dx > -\\ V N \\ L , \\G N \\ L2 > -C \\G N \\ L2 (\\G N \\ L2 + U N \ 



i" WL2 (Tin)) ■ 

(4.14) 
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The remaining part of (|4.13[) is also not very difficult. With the first integral on the r.h.s we have 



J n a N (d Xl a N + U -Va N )dx = \ J n (d Xl \a N \ 2 + U ■ V\a N \ 2 ) dx 
-\L \a N \ 2 dS+\L \a N \ 2 dS + U r U-n\a N \ 2 dS-U a \a N \ 2 dwUdx 
— ~ \ Ir a in dS + \ J r U ■ n\a in \ 2 dS 



U d ^ U hJW N \\l 2 



(4.15) 



>-c\W in \\l-E\\a N \\l 2 >-c-E(\\G\\i 

Hence we have 



\rr N \\ 2 



in\\L 2 



/ 3 > -C(||G||| 2 + \\Vin\\L 2 a\ J 



)-EU 



N\\2 



Combining ( |4. 1 l| ),( f47T2| ) and ( |4.16| ) we conclude 

(p N (t N ),n>c[u Nu2 



N 



D 2 ] 



(4.16) 



(4.17) 



where /J = ||F|| L2(n) + ||G|| L 2 (n) + ||#|| L2(n) + ||(7 in || L2 ( rin ) +Yh=i ||-Bi|U 2 (r)- Thus there exists 
M = M(/jl, ft, D) such that (P N (£ N ),£ N ) > for U N \\ w i = M, and applying Lemma 
g we conclude that 3£ N * : P N (£ N *) = and \\£ N * 
u n _ ^n* is a so i ut i on to g6|. 



| w i < M. By the definition of P 



N 

5 

□ 



Now showing the existence of weak solution is straightforward. The result is 

Lemma 16. Assume that F,G,H G L 2 (Q), a in G L 2 (T in ), E>i G L 2 (r), % = 1,2. Let \\U\\ W 2 
be small enough and assume n,L,f satisfy the assumptions of Theorem^ Then there exists 
(u,a,rj) G V x L 00 (L 2 ) x W 2 , that is a weak solution to system (3.1). Moreover, the weak 



solution satisfies estimate (3.4) 



Proof. Let us set a = S(G — divw ) and rj = R(u ), where u is the solution to (4.6). 



Estimates d2.26b, d4.5b and d4.8b imply that \\u N \\ m + \\a N \\ Loo{L2) + \\v N \\m < C(D Hn ). Thus, 



at least for a chosen subsequence (denoted however in the same way) 



.iV 



u 



u 



in H\ 



a 



.N 



^* a in L 00 (L 2 ) and r] N — ^ rj in W 2 



for some {u,a,i]) G H 1 x Loo(L 2 ) x H 1 . Passing to the limit in (|4~2[)— (|4~4|) for (u N ,a N ,i] N ) 
we conclude that (u,a,r]) satisfies ( |4.2 )-(4.4), thus we have the weak solution. To show the 
boundary condition on the density we can rewrite the r.h.s of (13. lh 6 as 



a 



G'-divt 
1+u i ~ v u V) 

on r, : 



in V'(tt), 



(4.18) 



and, treating x\ as a "time" variable, adapt Di Perna-Lions theory of transport equation (Q) that 



implies the uniqueness of solution to (4.18 ) in the class L 00 (L 2 ). The proof is thus complete. □ 
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4.2 Strong solution 



With the estimate ( |3.1 1[ ) the only problem is to tackle the singularities at the junction of T with 
T in and Y out . To this end we reflect the weak solution in a way which preserves the boundary 
conditions and then apply the classical elliptic theory to the extended solution of the Lame sys- 
tem. The details are given in Lemma [7] To show the regularity of the temperature we can apply 
Lemma [6] since on Y in U T out the boundary condition on the density reduces to the Neumann 
condition. 



5 Bounds on the approximating sequence 



In this section we will show the bounds on the sequence (2.24). Due to the term u ■ Vcr in 
the continuity equation we are not able to show directly the convergence of this sequence in 



Wp x Wp 1 x Wp to the strong solution of (2.20). We can show however its boundedness in W£ 



Wp x Wp. Next, using this bound we derive the Cauchy condition in W\ x L OQ (L 2 ) x W\, and 
thus show the convergence in this space to some (u, a, r/). On the other hand, the boundedness 
implies weak convergence in Wp x Wp x Wp, and the limit must be (u, er, 77); hence the solution 
is strong. 

The following lemma gives the boundedness of (u n , a n , rf 1 ) in Wi 2 x W p l x Wp. 

Lemma 17. Let {(u n , a n , rf 1 )} be a sequence of solutions to (2.24) starting from (u°, a°,rj°) = 
([0,0,0], 0,0). Then 

h n \\w2 + h n \\m + \\t\\-m <M, (5.1) 



where M can be arbitrarily small provided that D defined in (1.14), quantities \\B. 



i 1 5 2, || (Jin 

Theorem\l\ 



l \\w^{v in ) and 



are small enough and a, L, k, fulfill the assumptions of 



Proof. Estimate ( |3. 1 1| ) for system (2.24) reads 



\u 



n+l| 



1 ||_n+l|| 

w p 2 + If II wj 



\V llw| < 



< C [\\F(u n ,<T n ,7i n )\\ Lp + \\G(u n ,<j n ,ri n )\\ wi + ||FK,a",r/ n )|| Lp 
+ Z)i=i ll^llw^-^fr) + Ikmllwp 1 ^)]- 



(5.2) 



Denoting A n = \\u n \\ W 2 + ||cr n || w i + ||?] n ||iy2, from (2.23) and (5.2) we get 

A n+l <C(A 2 n + Al + D ). (5.3) 

We aim at showing that A n thus A n is bounded by a constant that can be arbitrarily small provided 
that Aq and D are small enough. Indeed, let us fix < 5 < (Note that the constant C can 
be without loss of generality taken larger than 1.) Assume that CD < 5. Then (5.3) entails an 
implication A n < 25 =3- A n+ i < 25 and we can conclude that 

\\u n \\ W 2 + IK'H^i + \\r] n \\ W 2 < 25 Wn G N. (5.4) 



□ 
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The next lemma almost completes the proof of the Cauchy condition in W\ x L oc (L 2 ) x W\ 
for the iterating scheme. 



Lemma 18. Let the assumptions of Lemma 17 hold. Then we have 



u 



n+l 



„.m+l\\ I _ n - m + 1 \ 



< E(M) (\\u n 



n+l 



^m+lll 

V WW} 



1l m \\ 1 



i II _n ^-mH i ||„n ^."il I ^ 



where M is the constant from (5.1) and E{M) can be taken arbitrarily small. 



(5.5) 



Proof. Subtracting (2.24) m from (2.24 | n we arrive at 



9 Xl (u""*" 1 - u m+1 ) - S(V(w n+1 - u m+1 )) + Pl V{cx n+l - a m+1 ) 
+p 2 V(r/ n+1 - r] m+1 ) = F(u n , a n , r] n ) - F(u m , a m , ri m ), 

div {u n+1 - u m+1 ) + d Xl (a n+1 - a m+l ) + {u n + w ) • V(cx n+1 - a m+1 ) 
= G{u n , a n ) - G(u m , a m ) - {u n - u m ) ■ Va m+1 , 

rid Xl (r) n+1 - r] m+1 ) + r 2 div {u n+1 - u m+1 ) - nA(r] n+1 
= H (V\ a n , 7] n ) - H(u m , a m , r/ m ) 

S(V(u n+1 - u m+1 )) -Ti + a {u n+1 - u m+1 ) ■ Ti\ r 

= 0, 
0, 



0, 



ra- (u n+1 -M m+1 )| r 

(T n+1 -CT m+1 | r . = 

-L 7,77. 



dn 



+ L(r/ n+1 - ri 



Estimate ( |3.4[ ) applied to this system yields 



L.n+l _„m+l|i I n„ 

I" " 1 1 Wo 1 ~r ll* 7 



n+l 



(7 



m+1 1 



0. 



Loo(L 2 ) + \\V n+1 ~ V 



n m+l\ 



Wi < 



|F(M n ,a n ,r/ n ) -F(M m ,a m ,r7 m )||v/* + ||G(w n , a") - G(u m , a m )\\ L2 

L2 + \\H{u n ,cr n ,r] n )-H{u m ,cr m ,r] m )\\ L2 . 



+ \\(u n -u m ) ■ Va 



m+l I 



(5.6) 

(5.7) 
(5.8) 

(5.9) 



In order to derive (5.5 ) from the above inequality we have to examine the r.h.s. The differences 
in G and H are bounded in a straightforward way: using the imbedding Wp ^ L ra and Holder's 
inequality we derive 



\\G(u n ,a n ) -G(u m ,a m )\\ L2 + \\H(u n ,a n ,r] n ) - H(u m ,a m ,r] m )\\ L2 
< E(M) (|K - u m \\ w , + Ha" - a m \\ Loo{L2) + ||rf - v m \\w})- ' 



(5.10) 



The difference in F must be investigated more carefully. A direct calculation yields 

F(u n , cr n , rf) - F(u m } a m , r] m ) = Ff' m + F 2 n ' m , 

where in F™' m we include all the terms which do not contain V(<r n — a m ) as well as other 
terms containing the gradient of density or its difference. Direct calculation using the imbedding 
H 1 L 6 , Wp <^-> and Holder's inequality yields 



\F?' m \\v* < E(M) (\\u n - u m \\ w , + h n ~ ° m \\L x[L2) ). 



(5.11) 
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Next, 



■(d p n(a n + 1,00 + 0! + rf) - p x )V{a n - a r 



(d p ir(a n + l,9 + 9 1 + r] n ) - d p ir(a m + 1,9 + 9 1 + f))Va m (5.12) 

— : *2,1 + -^2,2 ■ 

hence we multiply by v E V and integrate. With the first 



We have to compute V* norm of F 2 
term we have 

| j n -{d p ix(o n + l,9_o + 9j + rf) - Pi)V(a n - a m ) ■ v dx 
< | j n {d p Tx(a n + 1,00 + 0! + rf) - p^a 71 - a m )divvdx 
+ \ j n Vd p 7c{a n + 1, 0o + 0i + v n y v{a n - a m ) dx\ 
< C{\\v n \\w} + h n \\w} + ll^o - T \\ w i + \\9xWwi) \\(a n - OIUMk. 

hence 

\\F^\\ V * < E(M)\\(* n - *™)\\ Loo{L2) . (5.13) 

We have used the fact that \ + | + - < 1 and the imbedding W\ <-» L 6 . The other term can 
be estimated similarly, only without integration by parts. Combining the estimates on F{ l,m ,F^ m 
and F 2 "2™ we conclude 

\\F(u n ,a n ,fi n )-F(u m ,a m ,r, m )\\v < i?(M)(||n"-n m || H , 2l + ||a"- ( T m || Loo{i2) + ||^-r / m || W2l ). 

(5.14) 

The only term that remains to estimate is (u n — u m ) ■ Va m . We emphasize that this is the term 
which makes it impossible to show the convergence in x directly. Namely, if we would 
like to apply the estimate ( 3.1 1\ to the system for the difference then we would have to estimate 
\\(u n — u m ) ■ Va m \\wi what can not be done as we do not have any knowledge about ||c m ||w2. 
Fortunately we only need the L 2 -norm of this awkward term, which can be bounded in a direct 
way as 



II K 



u m ) ■ Va m || L2 < 



u 



u m \\ La 



|Va r 



< G \\a \\wi \\u 



n m \\ , 



since q 



2p_ 
p-2 



< 6 for p < 3. We have thus completed the proof of (5.5 ). 



(5.15) 



□ 



Now, Lemma [IT] implies that the constant E(M) < 1 provided that the data is small enough 
and the starting point (u°, a , rf) = ([0, 0, 0], 0, 0). It completes the proof of the Cauchy condi- 
tion in H 1 x L 00 (L 2 ) x H 1 for the sequence (u n , a n , rf 1 ). 

hold for any starting point (u°, a°,r]°) small enough in Wp x 



Remark. Lemmas 



17 



and 



18 



x Wp, but we can start the iteration from ([0, 0, 0], 0, 0) without loss of generality. 



6 Proof of Theorem U 

In this section we prove our main result, Theorem [T] First we show existence of the solution 
passing to the limit with the sequence (u n , a n , rf 1 ) and next we show that this solution is unique 



in the class of solutions satisfying (1.15). 
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Existence of the solution. Since we have the Cauchy condition on the sequence (u n , a n , rf 1 ) 
only in the space W 2 x L 00 (L 2 ) X W 2 , first we have to show the convergence in the weak 
formulation of the problem ( |2.20 ). The sequence (u n , a n , rf) satisfies in particular the following 
weak formulation of (|2.24[): 



J n {v ■ d Xl u n+1 + S(Vu n+1 ) : Vf — p\o n+1 div v — p 2 r] n+1 div v} dx 
+ f r a{u n+1 ■ Ti ) (v ■ Ti) dS = f a F{u n , a n , rf) ■ v dx + J^Li Ir B ^ v ' r *) dS Vwel/ 

(6.1) 



a n+l u n -V(j)dx- / (T n+1 div u n+l dx 



(j)(G(u n ,a r 



divw 



dx 



I in 



V <j> e C°°(fi) : <P\r out = 0, where u n = [1 + (u n + u )\ {u 11 + u ) 2 , {u n + u )% and 

(1 + e 2 ) J n d Xl r] n+1 w + T oP2 f Q wdivu n+1 + nf n Vw- Vrf +1 
+ J r Lwr] n+1 dS = J n H(u n ,a n ,r] n )wdx 

Now using the convergence in W 2 l X L^L^ x W 2 combined with the bound (5.1 



<7,; n dS 

(6.2) 



(6.3) 



Wp x Wp we can pass to the limit in (6.1 ) — { 6.3b. The convergence of the l.h.s. of (6T 
obvious. Recalling the definition (2.21 for F(-), G(-) and H(-) we verify easily that 



) in W3 x 



)-(6.3) is 



G(u n ,a n ) ■ v dx ->■ / G(u,a)-vdx 
u Jq 



(6.4) 



and 



/ H(u n ,a n ,ri n ) ■ v dx -> / H(u,a,r)) ■ v dx, (6.5) 
Jn Jn 

and the only step in showing the convergence of F which requires more attention is to show that 



(d p7 r(a n + lA+#i+?f)-pi)Va n W:r -> / (<9 p tt((7+1, ^o+^i + r/) -j?i)V(7-y dx. (6.6) 
n Jn 

However, since o n — > a and rf 1 — > r\ in C(f2) as well as Va" — 1 Va in L p (f2), we easily verify 



that ( |6.6[ ) holds true. We conclude that (m, <t, r\) satisfies (6.1 )-(6.3 ). Now we need to show that 



this implies the strong formulation (2.20), what can be done in a standard way, just integrating 
by parts in the weak formulation. 
Now we set 

v = v + u + u , p = a + l, 8 = r] + 6 + 9 1 

and we see that (v, p, 9) solves We clearly have E(M) = E(D ) where D is defined in 

( |1.14[ ), hence estimate ( |1.15[ ) holds. 



Uniqueness. The uniqueness in the class of solutions satisfying (1.15) actually results di- 
rectly from the method of the proof, more precisely from the proof of ( |5.5[ ). Namely, we can 
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show uniqueness on the level of perturbations, i.e. uniqueness for (2.20). Consider two solutions 
with the same data, denote it by (u\, ci, 771) and {u 2 , c 2 , ry 2 ). Their difference then satisfy 

d xi (ui - u 2 ) - S(V«i - Vw 2 ) + piV(<Ji - a 2 ) 
+P2^(Vi ~ V2) = F(u l7 0-1,171) - F(u 2 , cr 2 , 7/2), 

div («! - u 2 ) + ^(o"! - a 2 ) + («i + w ) • V(cr 1 - er 2 ) 
= 0"i) - C7(w 2 , CX 2 ) - Ox - m 2 ) • Vct 2 , 

(1 + e 2 )(9 :Cl (r/i - 772) + r p 2 div («i - u 2 ) - kA(t)i - r/ 2 ) 
= H(ui, <ri,r}i) - i?(« 2 , o"2, 7/2) 



S(V(«i - 7i 2 )77) -Ti + a(Ui - U 2 ) ■ Tj \t 

n- (ut - u 2 )\ r = 0, 
o\ - o 2 \ Vm = 0, 



0. 



dn 



' "" +L(r}i-r} 2 ) 



0. 



Now recall that to show < \5.5} we applied only ( ]5 .6^ — ( f5T9| ) and ( |5.1[ ), hence from the above equa- 
tions we conclude 



(6.7) 



\\U1 - U 2 \\ W 1 + ||0"i - 0- 2 || Loo(L2 ) + ||77i - 7] 2 \\ W 1 

< E(M) (\\ Ul -u 2 \\ w i + Wax -a 2 \\ Loo (L 2 ) + \\rji ~ mWw*)- 
Provided the data are small enough we have E(M) < 1 and so 

IK - u 2 \\ w i + |K - a 2 \\ Loo(L2) + || 77! - 77 2 |K 2 i = 
which completes the proof of uniqueness and hence of Theorem [TJ 
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